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ABSTRACT 

The Mg II h&k doublet are two of the primary spectral lines observed by the Sun-pointing Interface 
Region Imaging Spectrograph (IRIS). These lines are tracers of the magnetic and thermal environment 
that spans from the photosphere to the upper chromosphere. We use a double gaussian model to fit the 
Mg II h profile for a full-Sun mosaic dataset taken 24-Aug-2014. We use the ensemble of high-quality 
profile fits to conduct a statistical study on the variability of the line profile as it relates the magnetic 
structure, dynamics, and center-to-limb viewing angle. 

The average internetwork profile contains a deeply reversed core and is weakly asymmetric at h2. In 
the internetwork, we find a strong correlation between h3 wavelength and profile asymmetry as well hi 
width and h2 width. The average reversal depth of the h3 core is inversely related to the magnetic 
field. Plage and sunspots exhibit many profiles which do not contain a reversal. These profiles also 
occur infrequently in the internetwork. We see indications of magnetically aligned structures in plage 
and network in statistics associated with the line core, but these structures are not clear or extended in 
the internetwork. The center-to-limb variations are compared with predictions of semi-empirical model 
atmospheres. We measure a pronounced limb darkening in the line core which is not predicted by the 
model. The aim of this work is to provide a comprehensive measurement baseline and preliminary analysis 
on the observed structure and formation of the Mg II profiles observed by IRIS. 


1. INTRODUCTION 

The magnetic structure and thermodynamics of the so¬ 
lar chromosphere are open problems in solar physics with 
broad implications for stellar atmospheres. The visible 
spectrum contains several lines (H<a, Ca II H&K) which 
have allowed us to probe the structure and dynamics of the 
chromosphere using ground based observatories. The re¬ 
cent l aunch of the Interface Il egion Imaging Spectrograph 
(IRIS, |De Pontieu et al.||2Q14 ) provides a new dataset ob- 


serving one of the most important radiators of the chro¬ 
mosphere, the Mg II h&k lines at 2803.5A and 2796.4A, 
respectively. While there have been numerous missions 
to observe the Mg II h&k lines the measurements have 
been sparse. The earliest mea surements were done using 
a rocket-borne spectrographs (Durand et al. 1949). The 


Lemaire & Skumanich ( 

1973J). The OSO-8 (Artzner et al. 

T«)7T; 

vided 

I and Skylab (Doscl 

[ the first orbital di 

iek& Feldmanp977) missions pro- 

atasets. The Solar Maximum Mis- 


sion made th e first polarization mea surements of the Mg 


line profiles (Henze & Stenfio 1987). The highest spec¬ 


tral resolution measurements of the profile prior to IRIS 


were made by the HRT S rocket (Morrill et al. 2001) and 
the RASOLBA balloon (|staath & Lemaire 1995]) spectro¬ 
graphs. Center-to-limb measu r ements were d iscuss ed by 
Gouttebroze fc Lemaire l (|1974|), |Bonnet| (|1981|), and |Mor- 
rill & Korendyke (2008]). 'These measurements have been 
used to constr uct model atmospheres fo r a variety of so¬ 


lar structures. Gouttebroze (1977 1989) modeled the ef¬ 


fect of velocity gradient s in the internet work c hromosphere 
on profile asymmetry. Lemaire et al.| (1981) synthesized 
profiles for plage. Umbral profiles, which were originally 
identified as unique l y single peaked, hav e been studi ed by 


man (1984). Limb observations in _Mg II also provided 
diagnostics of prominences (Vial et al.||1979). 

In addition to the data points provided by solar observa¬ 
tions, multiple missions have conducted studies of stellar 
emission in the Mg II lines. The shape of Mg II h&k varies 
significantly across stellar types: from pure absorption in 
Altair to single peaked em ission in e Eri (Blanco et al. 
1982 Basri fc Linskyp979 ). The differences between the 


shape of Mg 11 h&k are a strong diagnostic of stellar winds 
and shocks in extended chromospheres like that of a Ori 

Variability and activity cycles 
re Mg II h&k lines in other stars 


(Uitenbroek et al. 1998) 
have been detected in t 


( Dempsey et al.|[l9 96). While stellar chromospheres exist 
in some form across the cooler half of the Main Sequence, 
there is a complicated relationship with coronae: a chro¬ 
mosphere is a prerequisite to form a corona, but stars wit h 
a chromosp here do not need to have a corona (Linsky & 
Haisch||1979|). The chromosphere-corona link is an impor¬ 
tant one for solar-terrestrial studies. Solar irradiance is a 


driving factor in determin ing the ion population s of the 
Earth’s upper atmosphere (Solomon & Qian 2005). While 
Mg II h&k are unlikely to play a significant role in thermo- 


Kneer et al. (1981), Lites & Skumanich (1982), and Gur- 


Mg II h&k are unlikely to play a significant: 
spheric photoionization, Mg II irradiance has been shown 
to be a sup erior proxy to EUV irradiance ove r the F10.7 
radio index ( Lean et aI7||2Q09 Guo et aL||2007| . 

The complex profiles observed in the solar Mg II doublet 
can be attributed to the lin es’ complex formation. A sem¬ 
inal study was presented in Milkey & Mihalas (1974). The 
Mg II doublet are resonance lines for a highly abundant 
element. As such they are very optically thick at the line 
cores, which form at relatively low densities under NLTE 
conditions. The formation depth of the profile varies sig¬ 
nificantly with wavelength. The core is estimated to form 
between 1-3 Mm above the (rsoonm = 1 ) photosphere. At 
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Structure 

Profiles (% 2 <1.5) 

Internetwork 

3.8 (3.3)xl0 6 

Network 

2.1 (1.4) xlO 5 

Plage 

1.2 (0.72) xlO 5 

Sunspot 

3.8 (3.4) XlO 4 

Filament 

1.4 (1-2) XlO 5 


TABLE 1 

Number of profiles in each structural region of the 24-Aug-2014 
dataset. Profiles with a qualifying % 2 are listed in parentheses. 

±0.5A the line opacity is reduced such that the formation 
layer maps near t he temperature m i nimum, 500 km above 
the photosp here. Leenaarts et al.| (2013a|b) and Pereira 
et ah (201%. used a ra diative 1VIHD numerical simulation 


(Gudiksen et al.||2011) to forward model the emission of 
the Mg 11 lines. These authors describe in detail the for¬ 
mation mechanisms of the lines, and how line components 
can be transformed into practical diagnostics of the solar 
atmosphere. Our research is complementary to these pa¬ 
pers. In this paper, we delve into the profiles of Mg II h 
observed by IRIS. We fit the line profile with a double gaus- 
sian model. A post processing routine parses the best fit 
model into a variety of profile statistics. We analyze a full- 
Sun mosaic dataset that offers us several million profiles, 
allowing us to present an overview of how the line com¬ 
ponents vary over the structurally distinct regions of the 
Sun. In particular, we analyze the variation of Mg II h with 
magnetic structure, dynamics, and viewing angle. Section 
2 presents the observations and data reduction. Section 3 
presents methodology for fitting the profiles and deriving 
the profile statistics. Section 4 describes the types of pro¬ 
files observed and interpretation on the physical sources of 
variability. Section 5 is summarizes the results and dis¬ 
cusses what future steps need to be taken. 


2. DATA 

IRIS is designed to investigate the heating of the solar 
atmosphere by collecting spectra and images of the Sun in 
three ultraviolet passbands: 1332-1358 A, 1381-1407 A, 
and 2783-2834 A. IRIS obtains spectral information by 
passing light from the Sun through a slit and onto a grat¬ 
ing. The spatial range of the spectra are limited by the 
175 arcsec length of the slit. Temporal limitations are in¬ 
troduced by the exposure time needed to build up suffi¬ 
cient photon counts, and by rastering the slit across the 
solar disk. 

We have designed an observing sequence that uses the 
capabilities of IRIS to build up a spectral map of the entire 
solar disk over as short a time as practicable. To gener¬ 
ate a full-Sun observation with the limited fleld-of-view of 
IRIS, we take successive observations at different satellite 
pointings and build up a mosaic of the Sun. 184 different 
observations were needed to generate the full-Sun mosaic. 
Each individual observation consists of a 64-step raster 
with 2 arcsec steps and 2 s exposure time at each slit posi¬ 
tion. The spectra along the slit have been binned to resolu¬ 
tion of 0.66 arcsec. This gives an area of 128 x 175 arcsec 
that takes ^ 190 s to observe. All 184 positions takes 
^ 18 hours, ^ 10 hours of which is observing time with 
the remainder in repointing the spacecraft. 

The above observing sequence is currently run approxi¬ 
mately once per month when IRIS is not in eclipse season. 
For the purposes of this study, we have chosen one such ob¬ 
servation from 2014-08-24 12:16 UT-2014-08-25 05:40 UT. 
The data used in this paper are IRIS level 2 data prod¬ 
ucts. These data have been processed from the raw ob¬ 


servations to remove bad pixels resulting from dust on the 
detector; dark current and flat field corrections have been 
implemented; and geometric and wavelength calibrations 
(based on the rest wavelength of Ni I 2799.17A) applied. 
For a complete description of the calibration process, see 
the IRIS user guid^] The full-Sun mosaic was assembled 
by positioning each raster according to its associated space¬ 
craft pointing. There are two sections of the mosaic that 
were observed while the orbit of IRIS was affected by the 
South Atlantic Anomaly. Due to an increase in cosmic ray 
hits at these times, we choose to remove these data from 
our analysis. 

This observation contains spectral information for sev¬ 
eral wavelength windows within the IRIS NUV and FUV 
passbands. For the purposes of this study, however, we 
limit our analysis to the Mg II h line at 2803.5 A. We an¬ 
alyze a spectral window 3.4 A wide, centered at the Mg 
II h line, with resolution of ^52 mA. In future work we 
intend to expand this analysis to include other emission 
and absorption features of the IRIS mosaic datasets. 

In order to compare the IRIS data with the photospheric 
magnetic field, we use simultaneous observations from the 
Helioseismic and Magnetic Imager (HMI, Scherre r et al. 
2012|) on b oard the Solar Dynamics Observatory (|Pesnell 
et al.|[2012|). HMI provides the line-of-sight component of 
the magnetic field for the entire Earth-facing solar disk at a 
cadence of 720 s. We combine many of these observations 
to arrive at a full-Sun line-of-sight (LOS) magnetic field 
map that most closely matches the times of observation 
of the IRIS data. For each IRIS raster we select the HMI 
observation that was temporally closest, and then select 
the portion of the HMI observation that matches the IRIS 
field-of-view. By repeating this process for all 184 IRIS 
rasters we can construct an HMI image of the full-Sun. 

3. METHOD 

We have chosen to use a 9-parameter double gaussian 
model to fit a 3.4A wide window centered on Mg II h: 


/(A) = a ± 6 * |A — c|±d* (exp(— — 2 ^ ) 

~h* exp( - )) (1) 

where the units of [/, a, d] are DN, b is DN A -1 , and 
[c, /? 9 : A fc] are A. Empirically, we find that this is an ac¬ 
curate model for a majority of the observed profiles. The 
Mg II h line is optically thick and the formation height of 
the profile varies by hundreds of kilometers from hi to h3. 
To calculate a theoretical profile, a synthetic atmosphere 
must be used to derive the source function as a function of 
wavelength along the line of sight. Our fit model does not 
attempt to include these physics, however the the double 
gaussian model provides a easily parameterized model that 
is capable of producing both single and double peaked pro¬ 
files with incongruous widths in the wing and core. The 
average solar Mg II h&k profile is an emission line with 
a reversed (depressed) core, which is captured by the su¬ 
perposition of a wide positive amplitude gaussian and a 
narrow negative amplitude gaussian. The far wings of the 
profile are captured by the linear function. While this tech¬ 
nique is applicable to Mg II k as well that analysis is com¬ 
plicated by the presence of the Mn I line at 2795. 6A. 


1 http://iris.lmsal.com/itn26/ 
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Fig. 1.— Typical double peaked profiles in the internetwork. Grey error bars show the IRIS data. The best fit model is plotted in dashed 
black. The BFM parameters for Equation 1 are listed in the legend. 
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Fig. 2.— Profiles with high x 2; non-gaussian inner wings (a), flat core (b), non-gaussian asymmetry (c), bright plage (d), cosmic rays (e), 
faulty minimization (f). 
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We use t he MPFIT least squares minimization algorithm 
(Markwardt 2009) to derive a best fit model. We assume 
tne measurement errors are the linear combination of Pois¬ 
son no ise (based o n the 18 photon per DN estimate in 


De 


Po ntieu et al.|2014 ) and a constant readout noise of 3 DN. 


After a best fit model (BFM, in terms of y 2 ) is retrieved 
from the minimization routine, we apply an algorithm to 
determine the number and location of local extrema in 
the BFM. The average profile will have two maxima at 
the red and blue edge of the spectral window at (see Fig¬ 
ure 1, at 170 km/s and -160 km/s respectively), bounding 
three minima (hlv, h3, hlr) and two additional maxima 
(h2v, h2r). The BFM value at these spectral positions 
is recorded along with the position. One additional case 
is accepted as a viable model: three maxima-two minima 
(single peak with weak wings). While some profiles are 
found with negative value for b in Equation 1, these pro¬ 
files are the result of noisy spectrums or high total y 2 . 

Incorporated with the IRIS package available for So- 
larsoft is the iris_get_mg_features algorithm which may also 
be used to identify the h2 peaks and h3 minimum (but not 
hi minima) of the h-line. We have opted for a fitting-based 
approach to derive a comprehensive parametrization of the 
profile that uses all the available data. A low-y 2 fit repre¬ 
sents an accurate model, while a high-y 2 fit can be further 
examined to determine how the model has failed. In a 
test sample, our fitting method and iris_get_mg_features 
differed by less than a spectral pixel for Xh 3 for 85% of 
profiles. 

We created structural masks to aid in parsing the 
dataset. The internetwork and network subsets were iden¬ 
tified using the SDO/HMI LOS magnetic field map. The 
magnetic qualifier for internetwork is that less than 5% of 
the pixels (0.5” pixels in the native resolution) within a 5 
pixel radius could have an unsigned LOS strength greater 
than 60 G. For network, more than 15% of pixels need 
to exceed that magnetic threshold. Plage was identified 
as large contiguous regions of unsigned magnetic flux with 


high intensities at 2803.3A (far enough from into the wing 
to be unaffected by the highly variable core intensity). The 
distinction between large network concentrations and plage 
was made manually based on the unipolarity of the nearby 
area. Sunspots were identified as contiguous regions of high 
unsigned magnetic flux with low intensities (<40 DN) at 


2801.8A (the outer most extent of our spectral window and 
the closest to continuum intensity). Filaments were identi¬ 
fied using only a spectral qualifier. Filaments were selected 
as contiguous regions where the line core intensity drops 
below 38 DN. Table 1 identifies the number of profiles in 
each region. 

Figure 1 shows typical profiles observed by IRIS in the 
internetwork. The extrema of the fit are labeled in the 
standard method: the reversed core is h3, the emission 
peaks are h2, and minimum in the wing is hi, v and r re¬ 
fer to the violet and red side of the profile. Our models fits 
return reduced y 2 values 0.4 ,which imply that our errors 
are a liberal estimate of the noise in the spectra. Approx¬ 
imately 48% of the on-disk spectra have fits with lower y 2 
than these models. We use y 2 < 1.5 as our viable model 
threshold, which applies to 86% of our dataset. 

The profiles in Figure la and lb have much in common 
and simultaneously many differences. Both profiles are 
double peaked, but the two peaks are not identical in in¬ 
tensity. To measure this difference we use the asymmetry 


statistic: 


Rh = 


Ih2v Ifi2r 


*-h2v 


Rh 


h2r 


[Pereira et al. (2013]). Figure la and 
i 0.13 and 0.15 respectively. The hi 


which is also used in 

lb have Rh values of 0.13 and 0.15 respectively, 
width (Xhir — Xhiv) and h2 width (A^r ~ A h 2 v ) are also 
similar for these profiles. However the core of the profiles 
are strikingly different. Both profiles have a h3 minimum, 
but Figure la has a much deeper core depression. We 
measure the relative depth of the core using the depth 
statistic: 

D h = 1 - 2Ih3 


lh2v 


Rh 


h2r 


Figure la and lb have Dh of 0.58 and 0.23, respectively. 
The differences of these profiles is emblematic of the com¬ 
plexity we observe in high resolution IRIS Mg II spectra. 
The correlations found in Leenaarts et al. (2013b) suggest 
that the intensity at the h3 core is inversely related to the 
altitude of the r = 1 layer, implying that Figure lb has 
less plasma in the upper chromosphere (pushing the r = 1 
layer lower) relative to Figure la. The extra chromospheric 
mass in Figure la could be related to dynamic mass loading 
from shocks or to dense overlying magnetic field. Further 
analysis on the properties of the h3 core will be discussed 
using maps of the h3 intensity in various magnetic struc¬ 
tures. 

Approximately 14% of the profiles for which we have 
sufficient data to fit have low quality BFM and are disre¬ 
garded for statistical analysis. These points occur through¬ 
out the dataset. There are several distinct types of faulty 
fits which are displayed in Figure 2. The most common 
bad fit look similar to Figure 2a. We have calculated the 
contribution to y 2 from three regions: far wing (window 
edge to hi), wing (from hi to h2), and core (h2 to h3). 
The error ratio for Figure 2a is approximately 2:7:6 (far 
wing:wing:core, respectively). The wing dominates the er¬ 
ror as the positive gaussian profile does not accurately cap¬ 
ture the changing slope of the profile between hi and h2. 
The model attempts to split the difference so the profile 
is too wide at hi and too narrow at h2. Of the 5xl0 5 
bad fits, 46% are dominated by wing errors. Figure 2b 
shows the second most common fitting problem: irregular 
cores. Many h3 cores tend to be flatter than the gaussian 
model. This is particularly true for low -/x profiles and in¬ 
side filaments. Approximately, 29% of the bad fit sample 
is dominated by core errors. 

In addition to flat cores, the double gaussian model is 
fundamentally limited in its ability to capture strong h2 
asymmetries. Figure 2c shows one such profile, where the 
model fits well the blue peak well but poorly the red peak. 
We find that 1% of the bad fits have a good fit for one peak 
and very poor fit for the other. That percentage is likely 
higher in actuality but is difficult to statistically identify. 

Some plage regions contain many bad fits. Indeed, ap¬ 
proximately 45% of the brightest plage regions (max[J] > 
800 DN) are bad fits. Figure 2d presents a typical bad fit 
plage profile. While the model captures the overall shape 
of the profile quite well, the small size of the error bars 
(primarily based on photon noise) magnifies the inadequa¬ 
cies of the double gaussian model. The bad plage profiles 
account for another 2% of the bad fits in the dataset. 

Figure 2e illustrates a cosmic ray hit data. We estimate 
that another 2% of bad fits are the result of anomalous 
spikes in the spectra. Figure 2f illustrates a fit where the 
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Fig. 3.— Overview of structure on 24-Aug-2014. Greyscale is the intensity at h3 where a successful fit was retrieved (left). Red pixels are 
single peak profiles. Green pixels are regions where a successful fit was not achieved. Blue pixels are pixels with insufficient data for fitting. 
The light blue boxes mark the regions displayed in Figure 4 and 5. SDO/HMI LOS magnetic field map (right). Filaments are bordered by 
yellow contours and plage is bordered by magenta. 


minimization routine honed in on a single peaked solution, 
where a double peak solution could also work. To eliminate 
this category of profiles, models that are single peaked but 
overestimate the intensity at the line center are flagged and 
the minimization routine is reinitialized with an narrower 
positive gaussian and deeper core. If the second BFM is 
still single peaked and overestimates the core intensity, the 
model is flagged. These profiles account for another 3% of 
bad fits. 


4. RESULTS 

4.1. Variation of Profiles based on Magnetic Structure 

Figure 3 shows the large scale structure of the Sun in 
h3 intensity and the magnetic field (zooming into the elec¬ 
tronic version of Figure 3 is suggested). The internetwork 
is largely faint, although there is significant small scale 
structure. Bright emission surrounds magnetic concentra¬ 
tions, network. Plage regions are very bright in h3. We 
do not expect the small flux tube effects of Spruit (1976) 
to apply at the estimated formation height of h3 (« 1 — 3 
Mm). Rather, the strong emission surrounding ma gnetic 
concentratio ns is related to stronger heating (Withbroe & 
Noyes|1977|) in the chromo sphere which has been w ell doc¬ 
umented in Ca II H&K (Skumanich et al. |1975 ). The 
dimmest emission occurs in filaments. 'The filaments we 
identify in h3 are similar to those visible in BBSO Ha 
fDenker et al.|[l999). Filaments are clear indicators of un¬ 
usual temperature stratifications in the solar a tmosphere 
(T < 10 4 K plasm a at altitudes of tens of Mm, Molowny- 


Horas et al. 


1999). The low temperature filament is likely 


scattering photons out of the LOS, thus reducing the in¬ 
tensity in a narrow Doppler broadened band surrounding 
the rest frame line center. 

Figure 4 and 5 display two enlarged regions from Fig¬ 
ure 3, Region 1 being quiet sun and Region 2 being plage. 
We have plotted a variety of statistics in these plots which 
measure the components of the profiles. The emission at 


hlv generally originates in the upper photosphere. Based 
on inspection of this and additional higher spatial resolu¬ 
tion IRIS datasets, we determine that the structures are 
a mix of reverse granulation and grains. At h2v, much of 
the small granular structure disappears and a diffuse halo 
of bright structures overlay and surround magnetic con¬ 
centrations. The brightest plage is almost twice as bright 
as the brightest network, and the maximum h2v emission 
occurs near the center of the magnetic flux concentration. 
Single peaked profiles are relativity common (^10%) in 
plage, as depicted by the contiguous red regions in Figure 
4. Figure 6a shows an example of a single peaked plage pro¬ 
file. While our algorithm separates between single peaked 
and double peaked profiles to aid in statistical analysis, 
there is in fact a smooth continuum of weakly separated, 
weakly reversed profiles that bridge those two profile cate¬ 
gories (Figure 6c shows an example). In the semi-empirical 
models of the solar atmosphere the Mg II h source function 
above the temperature minimum has a local maximum at 
an alt itude of 1.2 Mm (see Figure 7 in L 

a n 


Leenaarts et al. 
eight of 2 Mm. 


2013a), while the line core forms near 
4b produce a single peaked profile the atmosphere must 
be structured such that source function does not vary sig¬ 
nificantly between r = 1 at h2 and h3. At high densities 
the source function closely adheres to the Planck function, 
which rises through the chromosphere. These profiles are 
consistent with the model that hot high density loops are 
rooted in plage. Most sunspot profiles are single peaked, 
and an example is displayed in Figure 6e. Single pea ked 
umbral profile s we re previously reporte d in |Lites fc Sku- 


manich (1982) and Morrill et al. (120011). Sunspot profiles 


are half as bright at h2 and 20% narrower at hi than plage. 
The far wing intensity (beyond hi) is lower in sunspots 
than plage or internetwork. Single peak profiles also oc¬ 
cur in internetwork regions albeit at much lower frequency 
(^1%). Figure 6b shows a typical internetwork single peak 
profile. Although this profile contains a single maximum, 
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Fig. 4.— Regions extracted from Figure 3 displaying network and internetwork structure (left column) and plage (right column). The color 
coding is identical to Figure 2. The color table scaling is exponential with 7 = 0.25. The colar table bounds are listed in brackets. 


it still requires a substantial contribution from the negative 
gaussian amplitude (h in Equation 1) to achieve a low y 2 - 
Figure 6d and 6f show profiles that occur at low frequencies 
in the internetwork (< 0.0001%). These profiles represent 
some of the extreme variations from the mean profile shape 
that we observe in this large dataset. The profile of Fig¬ 
ure 6d could be produced by extreme upflows in the upper 
chromosphere which completely mask the photons in the 
nominal position of the h2v peak and wing. Figure 6f illus¬ 
trates Mg II h in absorption, which requires that the source 


function monotonically decrease from the photosphere pre¬ 
cluding the mid-chromospheric temperature rise. 

The intensity structure at h3 varies signifi c antly i n the 
quiet Sun. As discussed in 


Leenaarts et al. 


(2013b}, the 

Bifrost model predicts that there is a large range of for¬ 
mation heights for h3 (see Figure 4 in that paper). Semi- 
empirical models also predict higher column masses at a 
given altitude in the enha nced network versus faint inter¬ 
network (|Fontenla et al.||2009|), which alters the formation 
height for h3. The formation region of h3 in the network 
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Fig. 5.— Same as Figure 4. The color table scaling is linear 
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and in the internetwork likely varies significantly depend¬ 
ing on the magnetic topology of the upper chromosphere. 
Near the edge of network and extending into internetwork 
regions, coherent dark loop structures are visible. These 
are fibrils and there is a very strong contrast between the 
h3 intensity in and neighboring the fibril. Fibrils emanate 
from a fraction of network boundaries. In internetwork, 
it is more difficult to discern large coherent structures. It 
is likely that the chromospheric opacity in these regions is 
lower than in fibrils, and we are able to observe lower in 
the atmosphere. The structure at these lower heights is 


dominated by shocks so we do not observe large coherent 
horizontal structures. 

The h2v intensity and the hi width show similar struc¬ 
tures, both exhibit pronounced halos around flux concen¬ 
trations. If we consider an atmosphere where the Mg II h 
profile can be well described by a the linear+single gaus- 
sian model (Equation 1 with h = 0), then we could expect 
that there would be a strong correlation between the peak 
intensity and the width of the profile because we are mea¬ 
suring the width based on the inflection point at the con¬ 
tinuum and not the gaussian parameter g from Equation 
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Fig. 6 .— Single peaked profile in plage (a), single peaked profile in internetwork (b), weakly reversed plage profile (c), highly asymmetric 
single peaked profile (d), sunspot profile (e), reversed absorption profile (f). 


1. In the solar example where the profiles generally are not 
well described as a gaussian emission line, we would expect 
this effect to be reduced. In the internetwork (Solar-X be¬ 
tween -150 and -50 and Solar-Y between -500 and -440) , 
we see regions of broad hi width but weak h2v emission 
which provide a counter exam ple. The h2 width (also re- 
ferred to as peak sep aration in |Leenaarts et ah1|2013b and 


Pereira et al. 2013) shows significantly more small-scale 
structure. The h2 width can be correlated with the range 
of vertical velocities pr esent in the upper chromosphere 
(Leenaarts et al.||2013b). The internetwork h2 widths are 
not much narrower than in network or plage, despite a large 
difference in Ih 2 v In plage and network, collimated and 
linear structures are visible hinting at magnetically aligned 
features. The internetwork is likely traversed by complex 
web of canopy fields connecting ma ny dispersed and mixed - 
strength magnetic concentrations (|Schrijver & Title|2003). 
Flows along topologically distinct loops that cross through 
the LOS would not be correlated and might exhibit large 
velocity gradients that produce wide h2 separations. Con¬ 
centrations of low h2 width can be seen surrounding single 
peaked profiles in internetwork, network, and plage. As 
mentioned in the discussion of Figure 6, these profiles fit 


the gradual transition from weakly reversed double peak 
profiles to single peaked profiles. 

The Doppler shift of h3, A^, rarely exceeds 11 km/s. 
Near the edge of plage, fibrils show coherent velocity struc¬ 
ture. In the network, it is difficult to discern any fibrils in 
Aft 3 . Network actually exhibits a greater variability in h3 
velocity than plage or internetwork. The h3 velocity is 
highly correlated with the vertical velocity at the r = 1 
layer. The plage flow structures are likely related to the 
periodic flows and tran sition region oscillations reported in 


De Pontieu et al. (2003). 


4.2. Dynamic Profile Variations 

There are profound differences in the shape of the Mg 
II h profile and the morphology of spatial structure based 
on the magnetic environment. We also know the chromo¬ 
sphere is buffeted by waves driven by the turbulent velocity 
field of the photosphere. Many wave modes impart a di¬ 
rect modulation of spectral profiles based on the density 
and velocity perturbation. In addition, these waves may 
add to magnetic stressing or heating, which in turn varies 
the thermal structure of the chromosphere on small scales. 
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Fig. 7. — Example low -% 2 profiles from near disk center. On the y-axis are 100 stacked profiles sorted by the intensity at 2803.5A (v=0 
km/s). Each x-column represents profiles from regions defined by HMI LOS field strength. The spectral width of each column is ±55 km/s. 
The color table is exponential with 7 = 0.25 and represents the spectral intensity. Each profile with a well-defined h3 minimum has it marked 
with a blue bar. 


Based on these effects, we observe a high degree of variabil¬ 
ity of profiles within the magnetically classified regions. 

Figure 7 displays the mixture we see of magnetic and 
non-magnetic profile variability (zooming with electronic 
version of Figure 7 is suggested). Each column is sorted 
from lowest to highest intensity at the rest wavelength line 
core. Overall, there is a general trend that the brightest 
profiles occur in stronger magnetic field, which is consis¬ 
tent with the model that chromospheric heating may be 
induced or at least aided by magnetic fields. The bright¬ 
est profiles also tend to be single peaked. The intensity 
tends to be higher at h2v than h2r. This is an indica¬ 
tion that downflows predominate the chromosphere, which 
is well explained for shoc ks propagating through a non- 
magnetized atmosphere in Carlsson & Stein (1997). This 
effect has also previously been measured in opt i cally thin 
chromospheric emission lines by Peter & Judge (1999). 

The comparison between FigureH andTFIgure 7 is inter- 
esting. When structural information is included with the 
data (i.e. presenting in 2D maps), it is relatively straight¬ 
forward to categorize how magnetic profiles differs from 
non-magnetic profiles. By reducing the dimensionality of 
the dataset (i.e. presenting the profiles categorized only by 
magnetic field strength), the inherent variability of the pro¬ 
files becomes more obvious. Nearly all the types of profiles 
that we observe in the internetwork also occur in strong 


magnetic field regions, although the fractional occurrence 
rate may change. Table 2 presents examples of profiles 
with particular characteristics that are visible in Figure 7. 
The examples are selected to span both the magnetic field 
strength as well as profile core intensity. 

The variability depicted in Figure 7 is attributable to dy¬ 
namics as waves and shocks modulate the intensities and 
shape of the profiles. These processes should have a highly 
patterned and repeated spectral signature, similar to the 
Ha sawt ooth pattern attributed to passage of a shock in 
Hansteen et al. (2006). While our dataset is ill suited for 
determining those patterns, the temporal signature of the 
dynamics should be engrained in the distribution of the 
profiles. In Section 4.3, we present measurements on the 
mean and variance of a number of profile statistics. In 
addition to averaging the effect of dynamics into a single 
archetypal profile, it also important to consider the coor¬ 
dinated effects these dynamics have on the shape of the 
profile. To extract these relationships from our dataset, 
we calculate a cross correlation matrix of the profile com¬ 
ponents, presented in Table 3. The correlation coefficients 
allow us to measure if the variability of a particular statis¬ 
tic occurs synchronously with any other statistic. Here the 
cross correlation of statistics S a and S b is defined as: 


Cab — ^ 


(S? - Sa)(Si j ~ S b ) 




^a^b 
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Profile Characteristic 
Unshifted Single Peak 
Shifted Single Peak 
Wide hi 
Narrow hi 
Wide h2 
Narrow h2 

Asymmetric (I h 2v > Ih2r 
Asymmetric (I h2r < Ih3r 
Shifted h3 


Examples (Column and Rank) 
IN-86, 120G-99, 360G-8, 360G-99 
IN-15, 360G-30 

IN-28, IN-95, 200G-33, 320G-89 
IN-3, 240G-19, 240G-82 
IN-6, 200G-63, 240G-65 
80G-28, 80G-97, 160G-25, 360G-96 
) IN-43, 120G-28, 200G-86, 360G-75 

) 80G-55, 320G-14, 320G-72 

IN-73, 80G-55, 120G-34, 360G-77 


TABLE 2 

Notables profiles in Figure 7. Coordinates label the magnetic field 
and rank. 



Int. hlv 

Int. h2v 

Depth h3 

Asym. h2 

Width hi 

Width h2 

CO 

> 

Int. hlv 

1.00 

0.33 

0.11 

0.05 

0.22 

0.31 

0.06 

Int. h2v 

0.33 

1.00 

0.18 

0.43 

0.41 

0.10 

0.36 

Depth h3 

O.ll 

0.18 

1.00 

0.18 

-0.10 

0.46 

0.00 

Asym. h2 

0.05 

0-43 

0.18 

1.00 

0.00 

0.01 

0.80 

Width hi 

0.22 

0.41 

-0.10 

0.00 

1.00 

0.64 

0.04 

Width h2 

0.31 

0.10 

0.46 

0.01 

0.64 

1.00 

-0.03 

Vel. h3 

0.06 

0.36 

0.00 

0.80 

0.04 

-0.03 

1.00 


where S% v where represents statistic n at spatial position 
xy , S n is the mean and cr n is the variance of statistic n. 
C a b varies between-1 (anti-correlation) and 1 (correlation). 
We have limited our profile subset to high-quality double 
peaked fits of near disk center internetwork so as to limit 
the scope of variability to fine structure and dynamics and 
not viewing angle or magnetic structure. The variables 
which are most well correlated are the h3 wavelength and 
and h2 asymmetry. Because the far wings of the pro¬ 
file extend deeper into the atmosphere and do not shift 
in wavelength dramatically due to the low sound speed, 
this correlation is expected based on our double gaussian 
model. The more misaligned the negative gaussian is from 
the positive gaussian, the larger the h2 asymmetry. Given 
the success rate of our fitting algorithm and the quality of 
the fits, we believe this effect is robust against the bias 
of the fit model. The hi width and the h2 width are 
also st rongly correlated. Based on a scaling law argument, 


Ayres] (1979) suggests that the ratio of hi and h2 is deter¬ 
mined by the vertical extent of the chromosphere and rate 
of chromospheric heating. We expect that deeper profiles 
(high Dh) would produce wider h2 separation and thus 
have a strong correlation. The relatively high correlation 
between h2 asymmetry and h2v intensity is an artifact. 
By comparing asymmetry with h2r intensity and h-line 
radiance, we find brighter profiles are generally more sym¬ 
metric. The preponderance of positive asymmetry profiles 
biases the sample. 

These correlations are best used in conjunction with 
rapid cadence raster scans that allow us to extract both 
spatial and temporal information on the evolution of pro¬ 
files across well resolved spatial structures. In particular, 
we need to further develop studies which link the variable 
chromospheric profiles with the underlying photospheric 
drivers. 


4.3. Center-to-Limb Variations 

The viewing angle changes the observed intensities and 
shape of the profile because the column densities and veloc¬ 
ity projections are dependent on the line of sight. Figure 
8 and 9 illustrate the how the viewing angle affects the 
profiles in our three well populated magnetic subsets using 
joint probability distribution plots. The internetwork is by 
far the largest subset and is evenly distributed in y. The 
center-to-limb trends (slopes) are negative in hi, h2, and 
h3 intensity (similar to optical continuum) and positive in 
hi and h2 widths. The h2 and h3 intensity distribution 
are more dispersed at center than limb. This may be re¬ 
vealing that atmospheric thermal anisotropies decrease as 
a function of altitude. As constructed, there is no way to 
account for the spread of these statistics (at a given an¬ 
gle) using the semi-empirical models. The asymmetry and 


TABLE 3 

Cross correlation coefficients for components of internetwork profiles. 
The ensemble of profiles only uses models with % 2 < 1.5 and p, > 0.75. 
Abbreviations: Int.-intensity, Asym.-asymmetry, Vel.-velocity. 

depth statistics are unaffected by varying y. 

Given the clustered distribution of magnetic flux, 
the distribution of plage and network fields are less uni¬ 
form statistical samples. The magnetic field strength in 
plage varies significantly from region to region and plays a 
role in the scatter of the Ih 2 v and Ih 3 distributions. Plage 
and network are brighter than internetwork and wider at 
hi. For nearly all statistics, the distributions are more 
dispersed than internetwork. The hi widths are well cor¬ 
related with viewing angle. The h2 widths have a wider 
spread with an extended low-width tail. Plage and network 
are less asymmetric and have shallower cores compared to 
the internetwork. In order to measure radiance in the line, 
we have summed the data between hlv and hlr. In radi¬ 
ance near disk center, network is a factor of 3 brighter than 
internetwork, while plage is a factor of 8 brighter. Network 
is more dispersed in h3 velocity than internetwork or plage. 
Given the strong flows often observed in active region coro¬ 
nal spectra, it is surprising that network exhibits a higher 
p ercentage of h i gh-vel ocity profiles than plage. 


Avrett et al. (2013) discussed the center-to-limb varia¬ 
tion of the IVlg II k profiles based on the predictions of 
a series of semi-empirical models of the solar atmosphere. 
In that paper, Model B represents the mean internetwork 
while Model F represents the bright network. Table 4 com¬ 
pares the spectral intensities of the model and our derived 
profile statistics. While the semi-empirical models capture 
the variation of the emission peak (between internetwork 
and network and center-to-limb), it does not capture the 
variation in the core. The semi-empirical model predicts 
relatively little variation in the core intensity as a func¬ 
tion of /i, while we observe a measurable decline in the 
core at the limb. As discussed in section 4.1, if the forma¬ 
tion height of h3 is highly variable then the limb darkening 
effect can be explained by the increased sampling of high- 
lying and dark fibrils at low-/i. 

The semi-empirical models are constructed as a tool to 
incorporate a broad range of spectral information and con¬ 
vert it into a consistent model of the atmosphere. These 
static models provide us information on the structure of 
an averaged solar atmosphere, while the 3D MHD models 
provide information on the dynamics based on physics but 
limited by the computational power. By comparing the 
IRIS data with these models, we increase our understand¬ 
ing of this complex dataset and identify portions of the 
model which require adaption. 


5. CONCLUSIONS 
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Fig. 9.— Same as Figure 8. 


Statistic 

ALM 1 Model B 

ALM 1 Model F 

IRIS 2 IN 

IRIS 2 NW 

Peak Intensity, fi — 1.0 

210 

390 

170 

280 

Peak Intensity, fi — 0.1 

140 

270 

90 

- 

Core Intensity, ji = 1.0 

40 

130 

90 

200 

Core Intensity, /i = 0.1 

40 

120 

50 

- 

Depth, n = 1.0 

0.81 

0.67 

0.47 

0.29 

Depth, n = 0.1 

0.71 

0.56 

0.44 

- 


TABLE 4 

Comparison of center-to-limb variation from Avrett et al. (2013] ALM) for Mg II k and the measured quantities with IRIS for Mg II h in 
internetwork (IN and Model B) and network (NW and Model F). 1 ALM intensity units are W m -2 sr -1 A -1 . 2 IRIS intensity units of DN. 
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Most of our understanding of the chromosphere is based 
on ground based observations of Ca II and Ha. While 
datasets like that of SST/CRISP provide us the highest 
spatial resolution measurements of the chromosphere, they 
are limited in spatial FOV and the often brief periods of 
prime atmospheric seeing. The IRIS spacecraft has opened 
up a new window into the chromosphere, with the abil¬ 
ity to resolve sub-arcsec structure and the stability and 
atmosphere-free seeing to produce long-duration datasets. 
One of the primary purposes of the analysis of solar spec¬ 
tral lines is to diagnose the physical processes and the ther¬ 
modynamic conditions of the atmosphere. Prior to estab¬ 
lishing diagnostics, it is essential to identify the method¬ 
ology for measuring the necessary spectral features, and 
our works fits into this first step. The Mg II k line dif¬ 
fers slightly from Mg II h in formation height due to the 
difference in oscillator strength. It is important to extend 
this analysis to that line for additional diagnostic capabili¬ 
ties, however our method requires adaption. The Mn I line 
at 2975.5A can blend with the k-line for wide profiles. A 
well constrained fit model can likely be constructed with 
a spectral window that also includes the unblended Mn I 


2 999.1A line. This is a likel y dir ection for future an alysis. 
Leenaarts et al. ( 2013ajb ) and Pereira et al. (2013) have 


used an advanced numerical simulation of the solar atmo¬ 
sphere to forward model Mg II h&k emission and compare 
the derived spectral information which the physical proper¬ 
ties of the emitting plasma. Many of the statistics we have 
measured in this dataset are identical to those discussed by 
those authors. A direct comparison between the modeled 
profiles and observed profiles is difficult. While numerical 
models have proven powerful in describing specific physical 
processes (the formation of shocks by hi gh frequency ver¬ 
tical velocity fluctuations described the Carl sson fc St ein 


1997 for example), the models are not producing many of 
the - spectral features we identify in the observations. We 
know that the observed profiles are significantly broader 
than those forward modeled by Bifrost. Although the for¬ 
ward model incorporates magnetic fields, the limited scale 
and largely bipolar distribution of magnetic elements are 
no t direct analogues to a ll solar structures. As reported 


Leenaarts et al. ( 2012 ), the 3D radiation field (as op¬ 
posed to a ID vertical flux) can have a strong effect on the 
observed structure of the chromosphere. Ultimately, the 
distribution of Mg II ions, the relative populations of the 
ground and doublet states, and the incident radiation field 
along the line of sight are the determining factors for the 
shape of the emergent Mg II h profile. Both semi-empirical 
and MHD atmospheres provide us a means of calculating 
that distribution with strict limits on which physics can be 
included. 

The analysis of optically thick spectral lines present re¬ 


searchers with a double edge sword: the shape of the emer¬ 
gent profile encodes the thermal and hydrodynamic infor¬ 
mation of a significant vertical swath of the atmosphere, 
but the physical quantities we ultimately seek to measure 
are entangled. One potential diagnostic is a Mg II pro¬ 
file inversion, si milar to those used for optical spect ropo- 
larimetric data (Ruiz Cobo & del Toro Iniesta 1992). Of 
course inversions codes have largely operated in the regime 
where LTE is a reasonable approximation, which is not true 
for Mg II. The next step in this research is the integration 
of our derived quantities into our physical models (car¬ 
toon or otherwise) of the how the dynamic chromosphere 
is structured. 

In this paper, we have dissected a dataset containing 
Mg II h spectral profiles over the full solar disk. We de¬ 
composed the dataset into structural regions based on the 
magnetic field and quantified the differences between the 
regions based on a number of profile statistics. The inter¬ 
network is dim at h2 and deeply reversed. The network 
is brighter at hi, h2, and h3 than internetwork. It has 
enhanced hi width compared with internetwork. Plage is 
brighter and more variable at h 2 than network or internet¬ 
work. It has a small h 2 width and weak reversal. Single 
peaked profiles occur all over the disk, but most commonly 
in plage. We have created maps of many statistics for the 
line profiles. There are coherent but varied structures vis¬ 
ible in width, intensity, and velocity statistics. Fibrils are 
primarily visible in plage in X ^3 and in network in 1^ 3 . 
Both strong and weak field regions exhibit a great deal of 
profile variability which is the result of chromospheric dy¬ 
namics. We find strong correlations for the Doppler shift 
of the h3 core and the asymmetry in the h2 peaks and 
also between the widths of hi and h 2 . We calculated the 
center-to-limb variation of a number of statistics and com¬ 
pare them with the predictions of a semi-empirical model. 
We find that both h3 and h 2 exhibit limb darkening such 
that the core depth is flat as a function of /q while the 
semi-empirical model does not predict a strong limb dark¬ 
ening for h3. Our measurement of this effect is most robust 
in internetwork. This study presents the most comprehen¬ 
sive observational study on solar Mg II spectra to date, 
but it is clear that there are many facets on the formation 
and variability of this spectral line that are still not fully 
understood. It is essential to keep developing atmospheric 
models which can accommodate the great deal of dynami¬ 
cal and structural variability we observe in high resolution 
observations of the chromosphere. The photospheric ve¬ 
locity field strongly modulates the plasma conditions and 
radiation field in the lower chromosphere, and the upper 
chromosphere is highly structured by the magnetic field in 
the /3 < 1 regime. Both of these boundary conditions must 
be considered in our future modeling efforts. 
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